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Abstract 

We consider the problem of determining the unknown parameters of the Hamilto- 
nian of a network of spin ^ particles. In particular, we study experiments in which the 
system is driven by an externally applied electro-magnetic field and the expectation 
value of the total magnetization is measured. Under appropriate assumptions, we prove 
that, if it is possible to prepare the system in a known initial state, the above experi- 
ment allows to identify the parameters of the Hamiltonian. In the case where the initial 
state is itself an unknown parameter, we characterize all the pairs Hamiltonian-Initial 
State which give the same value of the magnetization for every form of the driving 
electro-magnetic field. The analysis is motivated by recent results on the isospectrality 
of Hamiltonians describing Magnetic Molecules. 



1 Introduction 



In recent years, chemists have developed methods to synthesize large organometallic molecules 
which contain a core of magnetic transition metal ions interacting via electronic superex- 
change interactions [3], [4], [5]. One of the main advantages of this technology is that it 
makes it possible to arrange the molecules in regular van der Waals crystals in which the 
magnetic interactions among different molecules are negligible. Thus, every cluster of this 
kind behaves like an assembly of identical and independent nanosize magnets, each corre- 
sponding to one molecule. For this reason, these novel systems are now deemed ideal to 
study fundamental questions concerning magnetism at the molecular level. In fact, even the 
simplest of these systems displays several new classical and quantum mechanical phenomena. 
One example is macroscopic quantum tunneling of magnetization [7] , [8] , a fascinating issue 
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of relevance to a variety of mesoscopic systems and nanostructures. This paper is a study on 
the determination of the parameters for Hamiltonians describing these magnetic molecules. 
The determination of exchange constants for these molecules has been traditionally obtained 
by measuring the temperature-dependent magnetic susceptibility or other thermodynamic 
properties of these compounds. This technique relies on the assumption that there is a one 
to one correspondence between the spectrum of the Hamiltonian of the system in a static 
magnetic field, and its thermodynamic properties as well as the numerical values of the 
model parameters. However, a recent article [13] has shown that quantum Heisenberg spin 
systems, such as magnetic molecules, may have different coupling parameters but the same 
energy spectrum and hence the same thermodynamic properties. This raises the question 
of whether a dynamic technique could be used to identify the parameters of a spin network. 
In particular, as in NMR or EPR experiments, one could try to let the system evolve under 
the action of a driving time varying electro-magnetic field, measure the total magnetization 
and, from the measured value, infer the value of the parameters of the system. The question 
of whether these experiments are able to identify the unknown parameters can be tackled at 
different levels according to how much we assume known about the system under investiga- 
tion; whether, for example, we assume the system prepared in a known initial state or the 
initial state itself a parameter to be identified. 

In this paper, we consider networks of spin | particles and present a positive answer to 
the above question. We prove that systems which give the same input-output behavior for 
any given state are the same. This shows that, if we can opportunely prepare an initial 
state we can use the above scheme to identify the parameters of the system exactly. If the 
initial state of the system is unknown we prove t hat there are only two possibilities (up to 
permutations of the spins) for two pairs Hamiltonian-Initial State to give the same input 
output behavior. They are either the same or the exchange constants have opposite signs 
and the initial states are related in a way we shall describe. So, in this case, the given 
experiments identify one of two possible systems giving the observed behavior. 

The paper is organized as follows. In the next section, we give the basic definitions and 
state the problem we want to solve in mathematical terms. In Section 3, we give the two 
results above described and prove the first one. This gives us the opportunity to elaborate 
on the control theoretic concepts of observability and controllability and their role in the 
parameter identification problem. The proof of the second result is much longer and it is 
presented in Sections 4 and 5. Section 6 presents some conclusions and a discussion of the 
results. 

2 Definitions and Statement of the Problem 

We consider an Heisenberg spin Hamiltonian of the form 

H(t) := i(A + B x u x (t) + B y u y (t) + B z u z (t)), (1) 
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with 

A '■= —i Y%<l,k,l=l Jkl(Ikx,lx + Iky,ly + Ikz,lz)i 

(2) 

B v := -i(T,k=ilkhv), for v = x,y, or z. 

For a network of n spin | particles, the matrix Ik ivi ,...,k r v r is the Kronecker product of n 
matrices equal to the 2x2 identity except in the kj — th (j = 1, r) position(s) occupied by 
the Pauli matrix cr„., i>j = x,y,z. Recall (see e.g. [12]) that the Pauli matrices are defined 
as 1 • 1 

a - := 2(l o)' ^ := 2~G 0*)' ^^(o -l)- (3) 

We denote by 72. the set of possible values for the density matrix, i.e. Hermitian, positive 
semidefinite matrices with trace one. Let p(t, u x , u y , u z , po), the density matrix solution of the 
Schrodinger equation corresponding to the controls u x (t), u y (t), u z (t) and initial condition 
Po- We assume that it is possible to observe the expectation value of the total magnetization 
in the x, y, and z direction, namely: 

M v (t) := M v (t; p Q , u x , u y , u z ) := Tr(S^ OT p(t, u x , u y , u z , p )), (4) 

where 

n 

s tot = j2l kvi for v = x,y, z. 

k=l 

We study the possibility of distinguishing the parameters by a single measurement of one of 
the above outputs. More precisely, we denote by E = S(n, Jkhlk) a model described by the 
equations (1) and (2), and by (S,p ) = Jki, 7fe), Po) a model given with fixed initial 

state po- Thus, for given control functions u x , u y , and u z , M v (t; p ,u x ,u y ,u z ), for v — x, y, 
and z, are the corresponding output functions. The parameters Jki and 7^ along with the 
number n of spins (and the value of the initial state po) characterize the model. The question 
of parameter identifiability through a single measurement of the magnetization can be posed 
by identifying the models (or set of parameters) that give the same input-output behavior. 
We have the following definition. 

Definition 2.1 Consider two models £ and E'. We mark with a prime ' all the symbols 
concerning system £'. 

• E and £' are equivalent and we write 

E ~ £' 

if and only if n — n' and for any given common initial condition po and control functions 

u x ,u y ,u z , we have 

M v (t;po,u x ,u y ,u z ) = M! u (t;po,u x ,u y ,u z ), for v = x, y, z. 
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• Two pairs model-initial state (X,Po) and (Yf,fJ ) are equivalent and we write 

(E,po)~(E', p' ), 
if for all control functions u x , u y , and u z , we have 

M v (t;p ,u x ,u y ,u z ) = M^(t;p' ,u x ,u y ,u z ), for v = x, y, z. 

Definition 2.2 A model is controllable if by varying opportunely the control functions 
u x ,u y ,u z , it is possible to drive the evolution operator from the identity to any unitary 
matrix. 

For general quantum systems controllability can be checked by verifying the so-called 
Lie Algebra Rank Condition [10] which means that the matrices A and B's characterizing 
the dynamics (cfr.(l)-(2)) generate the whole Lie Algebra su(n) (or u(n)) where n is the 
dimension of the underlying subspace (2 n for the case of networks of spin |'s). To a network 
of spin | one can associate a graph whose nodes represent the particles and an edge connects 
two nodes if and only if the corresponding exchange constant is different from zero. In the 
case of spin networks with different gyromagnetic ratios the system is controllable if and only 
if the graph associated to the network is connected, and sufficient controllability conditions 
can be given for the general case [2]. 

Definition 2.3 A model is observable if there are no two different states which give the 
same output for every set of control functions. 

Observability can be checked by verifying that the Observability Space V is equal to su(n) 
(n again is the dimension of the underlying Hilbert space) [6]. Considered the matrix that 
characterizes the output, Sj OT in our case, the observability space V is the vector space 
spanned by the matrices, 1 

ad k B n ad% 2 ■ ■ ■ ad%iS T v OT v = x, y, z, (5) 

where jk £ {0, 1,2,3} and B — A, B\ — B x , B 2 = B y , and B 3 = B z . The following fact 
holds true [6]. 

Lemma 2.4 Controllability implies observability. 

In parametric identification problems, it makes sense to restrict ourselves to observable 
systems since the unobservable dynamics does not contribute to the output which is our tool 
to identify the system. Moreover, we want to check that observable systems which give the 

1 ad R T := [R, [R, ...[R,T]]] where the Lie bracket is taken k times. 
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same input-output behavior, namely that are equivalents, have the same parameters. We 
can state the following two problems. 

Problem 1 Characterize the classes of observable equivalent spin models. 

Problem 2 Characterize the classes of observable equivalent pairs mo del- initial state. 

We shall see that the equivalence classes in Problem 1 consist of a single element. We shall 
solve Problem 2 restricting ourselves to networks that have different gyromagnetic ratios (or 
for which the spins can be selectively addressed) and controllable. 

3 Main Results 

In the following, we shall always denote by p(t) and p'(t) two trajectories corresponding to 
the same controls u x ,u y ,u z for the models E, £', respectively. The following Proposition 
whose proof we relegate to Appendix A, will be used used in the proof of both our main 
results. We notice here that the proof although presented for the case of spin Heisenberg 
systems can be adapted to any bilinear finite dimensional quantum control system. 

Proposition 3.1 Let (S,p ) and (E',p' ) be the two fixed pairs. Then, the following are 
equivalent: 

(a) (£,po)~(£',p[,), 

(b) For all control functions u x (t), u y (t), and u z (t), we have: 

Tr(Fp(t)) = Tr(F'p'(t)), (6) 
for all F G V and F' G V, with F' constructed as F changing all Bi in B\ (see (5)). 

We now state the first of our two main results. 

Theorem 1 Let E(n, J khlk ) := {A, B x ,B y ,B z } and E'(n,J' kl ,i k ) := {A 1 , B' x , B' y , B' z } be 
two equivalent models. Assume one of them is observable. Then A = A' , B x y z = B' . 

Proof. From the equivalence of the models and specializing (6) of 3.1 to a common initial 
condition, we obtain 

Tr(Fpo)=Tr(F'p ), (7) 

for every p G 1Z. Therefore F = F' . From the observability assumption F (and F') span 
all of su(2 n ) (n here is the number of spins, which is assumed to be the same). Since F is a 
generic element of V, we have [A, F] = [A 1 , F'\ = [A 1 , F] and therefore 

[A - A', F] = 0. (8) 
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Since F spans all of su(2 n ), A — A' must be zero. Analogously one can prove B x>y ^ z = B' x 



Notice that, although presented for the case of spin networks, Theorem 1 holds for any 
finite dimensional quantum system and essentially says that two observable models with the 
same input-output behavior, for a every state, must be equal. The assumption of observ- 
ability can be checked by checking controllability and applying Lemma 2.4. Conditions for 
controllability of spin networks are given in [2]. 

We now consider a more difficult problem since we assume to have much less knowledge 
of the model to be identified. We assume not to know its dimension nor the initial condition. 
We perform black-box type of experiments on two pairs model-initial state and we obtain 
the same results. We investigate what can be said about the two models. We assume that all 
the gyromagnetic ratios 7&, (7^) are different. This fact implies that the (mild) assumption 
that the graph associated to the spin network is connected is equivalent to controllability 
[2]. We shall assume this. Under this assumption, we can easily rule out the case in which 
the responses of the systems are both identically zero. In this case (and only in this case) 
the corresponding initial density matrices are scalar matrices and nothing more can be said 
about two equivalent models. So we will assume that the two initial states are not scalar 
matrices. 

Before stating the result, we need to introduce some more notation. We denote by I (X e ) 
the subspace of the Hermitian matrices of dimension 2 n generated by Kronecker products 
that contain an odd (even) number of Pauli matrices (and the rest 2x2 identity matrices). 
Moreover, if n is a permutation of the set {1, . . . ,n}, we denote by P n the matrix which 
transforms Kronecker products of n 2 x 2 matrices according to the permutation 7r (cfr. [9] 
pg. 260), namely for every n— pie of 2 x 2 matrices K ± , K n we have 



Theorem 2 Let (S,p ) = Jki, Ik), Po) and (£',p ) = (E'(n', J' kl , 7^), p' ) be two fixed 

models whose dynamics and output are given by equations (1), (2), and (4). Assume that 
both models are controllable, that all the 7^ and 7^ different from each other, and that p and 
p' are not scalar matrices. Then the following are equivalent: 



(b) n = n' and there exists a permutation n of the set {1, ... ,n} such that 
L 7fe = t4(*)> 

2. denoting by n] k = min{7r(/), vr(/c)} 7 and nf k = max{7r(/), 7r(fc)}, for 1 < I < k < n, 
then either: 



x,y,z- 

□ 



P^Kx ®K 2 ®---® K n )P n = K Al) ® K A2) <8) ... (8) K %{n) . 



(9) 



(a) (E, po) ~ (£', 




(10) 
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or 

( J lk = -J\ 2 VI < I < k < n, 

) lk Ktmrfm ~ ~ ' (11) 

[ pi = p[ and p 2 = -p' 2 ; 

where pi and p 2 (Vesp. p' x and p 2y ) are ine components of p (Vesp. P 7T p' P 7T ) inl , 
X e , respectively. 

Equations (10), (11) say that, up to a permutation of the spins, the exchange constants 
are all the same or all opposite. In one case the initial conditions are the same in the other 
case the components in X Q are the same while the components in X e are opposite. The next 
two sections are devoted to the proof of Theorem 2. We recall that we have two standing 
assumptions in all the treatment. The models we are dealing with have different 7's and are 
controllable. 



4 Preliminary Results 

In order to prove the implication (a) =^ (6) we shall need some properties of equivalent pairs. 
We present them in this section with all the proofs in Appendix B. The following proposition 
says that equivalent pairs (Model-Initial State) must have the same dimension. 

Proposition 4.1 Let (S,p ) and (£',Pq) be the two fixed models. If they are equivalent, 
then n = n' and there exists a permutation n of the set {1, . . . , n} such that: 

!■ 7fc = ln(k) for all A; e {1, ... , n}, 

2. Tr(I kv p(t)) = Tr{I^{k) v p'{t)) for all t > 0, all k e {1, . . . ,n}, all v G {x,y,z}, and all 
possible trajectories pit) of (S,p ) and corresponding p'{t) of (S',Pq). 

The proof of the next lemma is not presented in the Appendix since it is just a notational 
modification of the proof of Proposition 3.1. 

Lemma 4.2 Let W and W be two Hermitian matrices of dimensions 2 n and 2 n ', respec- 
tively. If, for every trajectory p(t) of (S,p ) and corresponding trajectory p'(t) of (£',p' ) we 
have 

Tr{Wp{t)) = Tr(W'p'(t)), (12) 

then for every F, F := adB h adB j2 • • ■ adB Jr W, and corresponding F', F' := adsr ads'. • • • 
ad B > W, with r > and ji, ...,j r G {0, 1, 2, 3}, we also have 

3r 

Tr{Fp(t))=Tr{F>p'(t)). (13) 

Lemma 4.3 Let W and W be two Hermitian matrices of dimensions 2 n and 2 n ', respec- 
tively. If, for every trajectory p(t) of (S,p ) and corresponding trajectory p'(t) of (S',p' ) we 
have 

Tr{Wp{t)) = Tr(W'p'(t)), (14) 
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then, for every K, K := ads, ads, • • • ads, Bj , and corresponding K', K' := ads', ads 1 . ■ 
• ■ ads'. Bj Q , with r > and jo, ■■■■,j r G {0, 1, 2, 3}, we also have 

Tr({W, K]p(t)) = Tr([W, K']p'(t)). (15) 

Lemma 4.4 Let (£,po) and (X',p ) be two fixed models. Assume that they are equivalent 
and let 7r be the permutation given by Proposition 4.1. If W and W are two given Hermitian 
matrices such 

Tr{Wp{t)) = Tr(W'p'(t)), (16) 
for every pair of corresponding trajectories pit) and p'(t), then it also holds 

Tr([W, I kv ]p{t)) = Tr([W, h {k)v \p> \t)) , V k G {1, . . . , n}. (17) 



5 Proof of Theorem 2 

Let (S,po) and (E',p ) be the two given equivalent models. Assume that both models are 
controllable and that all the 7^ and 7^ are different from each other. We already know that 
n — n' and 1 and 2 of Proposition 4.1 hold. To simplify the notations, we assume, without 
loss of generality, that we have performed a change of coordinates in the second model so 
that the permutation 7r of Proposition 4.1 is the identity. Thus we can write 

7fe = 7fc, Vfce{l,...,n}, (18) 

and 

Tr(I kv p(t))=Tr(I kv p'(t)) V k E {1, . . . , n}. (19) 
Equations (10) and (11) now read as: 

f J lk = J' lk VI < l < k < n, 

1 p'o = Po; 

or 

f J lk = - J' lk VI < l < k < n, 
\ pi= p[ and p 2 = -p' 2 ; 

where p\ and p 2 (resp. p[ and p' 2 ) are the components of po (resp. p' ) in I , I e , respectively. 
We shall need the following two lemmas whose proofs are presented in Appendix C. 

Lemma 5.1 Assume that for all t > 0, all possible trajectories pit) of (S,po) an d corre- 
sponding p'(t) of (E', po), for fixed values 1 < ki,...,k r < n, and fixed Vj G {x, y,z} we 
have: 

Tr (I klVu ..., krVr p(t)) = Tr (I kl 

Vl,...,k r V r P'(t)) , ( 22 ) 

Then: 

1. equation (22) holds for any possible choice of the values of Vj G {x,y,z}; 
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2. 



1 J-klVl^.^krVr 



p(t))=Tr([[iI h ,[iI- kv - k ,A] 



1 1 klVl,...,krVr 



P'(t) 



for every values 1 < I ^ k < n and every {vj ^ v ^} G {x, y, z}. 



(23) 



Lemma 5.2 Assume that for all t > 0, all possible trajectories p(i) of (S,p ) an d corre- 
sponding p'{t) of (S',Pq), for fixed values 1 < ki, . . . , k r < n, Vj e and for given 
constants a and a', we have: 



aTr (Ik lVl ,...,k r v r p(t)) = a'Tr (/fc 1 « 1 ,...,fc rWr p , (t)) , 



(24) 



Then 



1. For any pair of indices fc, I e {1, . . . , n} with e {/ci, /c r } and / ^ {/ci, fc r }, 

aJlsTr (lk lVl ,...,k r v r ,TvP(t)) = a'J L k i Tr (h lVl ,...,krv r ,ivP'(t)) , (25) 
for any value v G {x, y, z}. 

2. For any pair of indices k, I both in {k\, k r }, (for example k — k±, I — A^) then 

aJjsTr {h^Mvz^krVrpit)) = a'JifTr (/k 1 « 1) fc 3 ,* ) ... ) fc r w r p , (t)) • (26) 



5.1 (a) 



Fix any 1 < k\ < k 2 < n, then, by applying statement 1. of Lemma 5.2, i.e. equation (25) 
with k = k\, I = k 2 to equation (19) with k = k±, we have: 



J klk2 Tr (I klVl ,k 2 v 2 p( t )) = A x k 2 Tr {hrv^Mvd {*)) , v vi, v 2 G {x, y, z). 



(27) 



Now, to the previous equality, we apply statement 2. of Lemma 5.2, i.e. equation (26) with 
k — k\ and I — k 2 to get: 



Jt lk2 Tr (h lVl p(t)) = j' k \k 2 Tr (h^p'it)) 
which, by equation (19), implies: 



72 _ j'2 
J kik 2 — J kl k 2 - 



(28) 



Therefore the exchange constants are equal up to the sign. Now we prove, by the way of 
contradiction that they are either all equal or all opposite i.e. 



{ J kl = J' kl V 1 < k < I < n} or { J kl = -J' kl \/l<k<l<n.} 



(29) 



Assume, by contradiction, that (29) does not hold. By the controllability assumption and 
the results of [2] we know that the graph associated to the network is connected. From this 
fact it is not difficult to see that if (29) is false, there must exist 3 indices I, ki, and k 2 (here, 
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to simplify notations, we assume 1 < I < k\ < k 2 < n, the other cases can be treated using 
exactly the same arguments) such that: 

J iki = J'lkv and J ^2 = - J 'ik 2 - (3°) 

Using equation (27) we get: 

Ji kl Tr(I lvMvi p(t)) = J' lki Tr{I lvMvi p'{t)), 

for all v,Vi G {x,y,z} and all corresponding trajectories p(t) and p'(t). By applying to the 
previous equality statement 1. of Lemma 5.2, i.e. equation (25) with k = I and I — k 2 , we 
get: 

W (31) 

for all v,Vi,v 2 G {x,y,z} and all corresponding trajectories p(t) and p'(t). Now we apply to 
equation (31) statement 2. of Lemma 5.2, i. e. equation (26) with k — k\ and I — I to get: 

'(*))■ (32) 

On the other hand, we can apply to equation (31) again statement 2. of Lemma 5.2, i. e. 
equation (26) this time with k = k 2 and I — I to get: 

Jik 1 Jm 2 Tr {h 1 v 1 ,k 2 v 2 p{ t )) = ■J'ik 1 J ik 2 Tr ( I k 1 v 1 ,k 2 v 2 p'( t ))- (33) 

Since Jf kl Jik 2 = -J'ilAk 2 while J «fci J ffc 2 = J'lkAh ( b ^ ( 30 ))> equations (32) and (33) imply 
Tr(lk ivit k 2 v 2 p(t)) = for all v 1 ,v 2 G {x, y,z} and all trajectories p(t). This fact contradicts 
the controllability assumption, thus equation (29) holds. 

If Jki = J' kl for every pair 1 < k < I < n, from the observability of the model, we must 
have po — Po, thus equation (20) holds. 

On the other hand, if J' kl = —Jm for every pair 1 < k < I < n, we argue as follows. First, 
we prove, by induction on 1 < r < n that: 

Tr(I klvl: __. tkrVr p(t)) = (-l) r ~ 1 Tr(4 1 „ li ... iA . r „ r p'(t)), Vuj G {x, y, z}, 1 < k x < . . . < k r < n, 

(34) 

and for all corresponding trajectories p(t) and p'(t). 

For r = 1 the previous equation is equation (19), thus the result holds for r — 1. Assume 
that (34) holds for 1 < r < n, and consider an arbitrary set of indices 1 < k 2 < . . . < k r+1 . 
By the inductive assumption we have: 

Tr(I klVl _ krVr p(t)) = (-l) r - 1 Tr(I klVl _ krVr p'(t)). 

Since the graph associated to the network is connected, for each node kj, j = 1, . . . , r, there 
exists a path joining the node kj with the node k r+1 . Let j be the index for which this 
path is the shortest one, and denote by h,...,ld the intermediate nodes. By the way we 
have chosen j, it is easy to see that {k±, . . . , k r } fl {h, . . . , la} — 0. To fix notations, we may 
assume without loss of generality (being all the other cases the same) that 1 < k-i < . . . < 
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k r < l\ < . . . < la < fc r +i- By applying statement 1. of Lemma 5.2, equation (25) with k — j 
and I — li, since Jj h = —Jj tl , we have: 

Tr(I klVlj _ tkrVrthwi p(t)) = (-l)(-l) r-1 T , r(4 1 „ lv .. jfcr „ 7 . )JlU)1 ^'(i)), 

for any w\ G {x, y, z}. By applying again statement 1. of Lemma 5.2, equation (25) another 
d — 1-times with k = k and I = i — 1, . . . ,d— 1, and then another time with k — Id and 
I — k r+ i, we end up with: 

T r (Ik 1 v 1 ,...,krVr,hw 1 ,...,l d w d ,k r+1 v r+ iP(t)) — ( — 1) ( — 1) Tr(^Ik 1 v 1 ,...,k r v r ,l 1 w 1 ,...,l d w d ,k r+l v r+1 P 

'(*))■ 

(35) 

Now we apply to equation (35), statement 2. of Lemma 5.2, i.e. equation (26) d — 1-times 
with fc = l i+ i and Z = Zj, i — 1, . . . , d — 1, to get: 

Tr(Ik lvl ,...,k r V r ,l d W d ,k r + 1 Vr+lP{t)) = ( _ 1) ( _ 1) ( _ l) r ^ r (^fclDl,...,fcrfr,idW d ,fcr + l«r + lP '(*))■ 

Finally, by applying again statement 2. of Lemma 5.2, i.e. equation (26) with k = k r+ i and 
I = Id, we end up with: 

Tr{I k 

p'{t)), 

as desired. Thus equation (34) holds. 

Now, denoting by p\ and p 2 (resp. p[ and p' 2 ) the components of p (resp. p' ) in X , X e , 
by using (34), we have: 

Tr((/ O (p -P , o)) = 0, Tr((I e (p + p' )) = 0, (36) 

for all elements I G X G and 7 e G X e . Equation (36) implies that the components of p and 
p' Q in X coincide while the components in X e are opposite to each other. Thus equation (21) 
holds. 



5.2 (6) => (a) 

Let (S,Po) and (^'>Po) be two models, which are both controllable, with all the ^ k and 7^. 
different from each other. 

In the case where equation (20) holds (i.e. same model and same initial condition), 
obviously that: 

p'(t)=p(t), 

for all t > 0. Thus the two models are equivalent. 
Assume now that equation (21) holds. Thus, 

A' = -A, and B' v = B v V v G {x, y, z). 

and 

p 1 =p' 1 and p 2 = -pa- 
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We have: 

p{t) = [A + B x u x {t) + B y u y {t) + B z u z (t),p(t)\ , (37) 

while 

p'(t) = [-A + B x u x (t) + B y u y {t) + B z u z {t),p'{t)\ . (38) 
It is easily verified (cfr. Lemmas 5.1 and 5.2) that: 



[B,I ] G 1 , [A,I ]El e , VI El , 

[B,l e ]el e , [A,l e ]el Q , VJ e eZ e . 

Thus, we can write the differential equations for pi(t) and p 2 {t) as: 

Pl (t) = \B x u x (t) + ByUy(t) + B z u z (t), Pl (t)} + [A,p 2 (t)} 
p 2 (t) = [A, Pl (t)\ + [B x u x (t) + B y Uy(t) + B z u z (t),p 2 (t)\ 

and similarly the differential equation for p[(t) and p' 2 (t) as: 

= [S x « x (t) + B„ti,,(0 + /»!(*)] - [A,p' 2 (t)\ 

p' 2 (t) = -[A,p[(t)] + [B x u x (t) + B yU y(t) + S^(t)y 2 (t)] 



(39) 



(40) 



(41) 



Combining equations (40) and (41), we obtain a differential equation for pi(t) — p[(t) and 
for p 2 (t) + p' 2 (t). In particular, we have 

Pl (t)-p[(t) = [B x u x (t) + ByUy(t) + B z u z (t),p 1 (t)-p' 1 (t)] + [A,p 2 (t)+p' 2 (t)] 
P2(t)+p' 2 (t) = [A, Pl (t)- p[(t)} + [B x u x (t) + B y u y (t) + B z u z (t),p 2 (t)+p' 2 (t)} y > 

From equations (42) it follows that if pi(0) = p[(0) and p 2 (0) = —p' 2 (0) then p 1 (t) = p[(t) 
and p 2 (t) = —p' 2 (t), for every t and for every controls u x (t), u y (t), and u z (t). In particular, 
since Tr(Sj OT p(t)) = Tr(Sj pi(t)) for all v G {x, y,z} and p\(t) = Pi(t) the two models 
are equivalent. 



6 Conclusions 

In this paper, we have investigated methods of dynamic parameter identification for networks 
of spin | particles. We have shown that by driving the network with an appropriate electro- 
magnetic field and measuring the total magnetization in a given (arbitrary) direction it 
is possible to identify the parameters. Moreover, if the initial state is not known, it is 
possible to obtain combined information about the initial state and the parameter values. 
We have assumed that all the gyromagnetic ratios of the spins are different or that it is 
possible to address each spin separately. In the opposite case, where all the gyromagnetic 
ratios are the same, the unitary evolution X(t), solution of Schrodinger operator equation, 
has the form X(t) = e A *<3>(£) where $(£) depends only on the controls u x ,u y ,u z and A is 
defined in (1), (2). In this case, we have that A commutes with $ and Sj OT , v — x, y, z, and 
therefore the output Tr{S^ OT e At ^{t)p^*{t)e- At ) is equal to Tr(e~ At Sj OT e At ^(t)p ^*(t)) = 
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Tr(Sv OT §(t)po$*(t)). The output is therefore independent of A. This implies that it is not 
possible to identify the parameters in A by a reading of the total magnetization. 

In our approach, the system theoretic concepts of controllability and observability as 
well as previously known results on the controllability of spin networks have played an 
important role. This is usually the case in the theory of parameter identification and we 
believe this approach will be useful for other classes of quantum systems. Extensions of 
the results presented here are possible and will be object of further research. For example, 
the hypothesis of controllability of the models can be weakened. If a spin network is not 
controllable and has different gyromagnetic ratios the associated graph has several connected 
components. The dynamical Lie Algebra associated to the system is the direct sum of 
Lie Algebras isomorphic to s-u(2" J ) where rij is the number of nodes (spins) in the j—th 
component [2] . Another important research problem is the actual design of control algorithms 
for parameter identification for which the research presented here is a preliminary step. 
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Appendix A: Proof of Proposition 3.1 

Proof It is clear that (b) implies (a), since iSj OT G V (iS'J OT G V) for v = x, y, z. We will 
prove the converse implication by induction on the depth s = J2i=i ki in (5) of the matrix 
F G V. If s = then F = iSj OT for v = x, y or z, thus equation (6) holds by definition of 
equivalence. Assume that equation (6) holds for matrices in V of depth < s and let F G V 
with depth equal s + 1. Then 

F = ad B -G, with G G V, 

and the depth of G is equal to s. Assume, by contradiction, that there exist control functions 
u x (-), Uy(-), and u z {-), and t > such that equation (6) does not hold, i.e. 

Tr{Fp{i)) ± Tr{F'p\t)). (43) 

On the other hand, since by the inductive assumption, equation (6) holds for the matrix G, 
we have: 

jTr(G P (t)) = jTr(G'p'(t)), 

for alH > 0. This implies: 

Tr([G, B ]p(t)) + Tr([G, B 1 ]p(t))u x (t) + Tr([G, B 2 ]p{t))u y {t) + 

Tr([G,B 3 ]p(t))u z (t) = Tr([G',B' }p > (t)) + Tr([G',B' 1 ]p'(t))u x (t) + 

+Tr([G', B' 2 ]p'(t))u y (t) + Tr([G', B' 3 ]p'(t))u z (t) (44) 

Define: 




u v (t) for t < t 
for t > i. 
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Then, clearly the trajectories p(t) and p'(t) corresponding to the two sets of controls u v (-) 
and are equal up to time t. Thus evaluating (44) at t — t , using controls u° v , we get: 

Tr([G, B ] P (t)) = Tr([G>, B' ]p'(t)), (45) 

which contradicts (43) if % — 0. Assume % ^ 0. First notice that, by repeating the same 
argument as above for a generic t > 0, the previous equality (45) must hold for all t > 0. To 
get a contradiction we use the control functions u l v given by 2 : 



Kit) 



u v (t) for t < t 
b\ v for t > t . 



Again the trajectories corresponding to the two set of controls u v (-) and u % v (-) are equal up 
to time t, thus evaluating (44) at t — t using controls u l v we get: 

Tr([G, B ]p(t)) + Tr([G, ^]p(t))^(t) = Tr([G>, B' ]p'(t)) + Tr([G>, B{]p' (t)) Ul (t) 

which, since the first terms are equal as observed before, contradicts (43), and ends the 
proof. □ 



Appendix B: Proofs of the preliminary results in section 
4 

Proof of Proposition 4.1 

Proof. It is not difficult to see that, for all I > 0, if F l — Ylk=ilk^kv £ V, then the 
corresponding matrix in V' is F' 1 = J2k=i l^kv By Proposition 3.1, it holds that: 

Tr(F l p(t)) = Tr{F' l p\t)). 

for all t > and all possible trajectories. Thus we have: 

n n' 

J2riTr(I kv p(t)) = £7?TV(W(*))- (46) 
fc=i fc=i 

Fix v and let aik(t) = Tr(Ik v p(t)) and cx' k (t) = Tr(Ik v p'(t)), then we rewrite equation (46) 
as: 

n n' 

EtWO-Et^W = (47) 
fe=i fe=i 

2 ^ix = = $3z = 1, <5i„ = otherwise. 
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The matrix M e ft(«+«')x («+«') given by; 



/ 



M = 



7i 

7? 



1 

In 
'n 



1 

Vi 

7? 



/n+n' 
7l 



1 \ 

7n' 

t3 



is a Vandermonde type of matrix. Notice that the coefficients 7&, k — 1, . . . ,n and also 7^, 
k — 1, . . . ,n', are all different. Moreover, the coefficients «*;(£) and a4(£) are not identically 
zero. In fact, if «*;(£) was identically zero we would have Tr(Ap(t)) = 0, for every A e sw(2 n ) 
(by the controllability assumption) which would imply p equal to a multiple of the identity 
matrix which we have excluded. Thus, from equation (47), we conclude that there exist two 
indices k and n(k) such that 

Ik = tC(S)- 

We can rewrite equation (47), as 



i («*(*) 



a 



vr(fc) 



+ E <Mt) - E tX(*) = 



(48) 



k=\,k+k 



k=l,k^n(k) 



Now we can repeat the same argument and, unless n = 1 or n' = 1, we will find two more 
indices j and 7r(j) whose corresponding values of 7j and 7^ are equal. We may assume 
without loss of generality that n' > n and repeat this procedure n-times. Thus we find a 
permutation n from the set {1, . . . , n} to the set {!,..., n'} and we rewrite equation (48) as: 



7r(fc) 



fc=l 



E 

fc ^ 7r(j) 

j = l,...,n 



tX(*) = 0. 



(49) 



Now, we can apply again the same argument, using the Vandermonde matrix N constructed 
with all the coefficients 7^ and the coefficients 7^ for those indices that are not in the image 
of 7r. Since the coefficients a' k (t) are not identically zero, we can conclude that all the 
coefficients j' k , for those indices that are not in the image of 71, must be zero. Thus, in 
particular n = n', the map n is a permutation and : 



7fc = 7^( fc ), <Xk(t) -a' <k) (t) = 0, 



Vt > 0, 



which concludes the proof. 



□ 
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Proof of Lemma 4.3 

Proof. We will prove the result by induction on the depth r of K and K 1 . If r = 0, the 
result follows from Lemma 4.2. Now assume that for every pair K and K' of depth r, and 
every pair of matrices W and W, (14) implies (15). From the inductive assumption we have 

Tr({W, K]p(t)) = Tr([W, K']p'(t)), (50) 

and 

Tr([W,BMt)) = Tr{[W\B' 3 \p'{t)). (51) 

Applying the inductive assumption with W (W) replaced by [W : K] and [W, Bj] ([W',K'] 
and [W,S^-]) we obtain 

Tr([[W,K],BMt)) = Tr([[W',K'],S>'(0), (52) 

TrmB^KW)) = Tr{[[W',B' 3 ],K']p'{t)). (53) 
Combining (52) and (53) using the Jacobi identity, we obtain 

Tr([W, [K,BMt)) = Tr([W, [K> \ Bflp? '(f)) , (54) 
which proves the lemma. □ 



Proof of Lemma 4.4 

Proof. For any Z > let Fi := YJk=\ llhv £ V; then its corresponding matrix in V' is 

F 



F{ '■= Sfc=i ll(k)Iir(k)v By applying Lemma 4.3, with K = F t and K' = F(, we obtain 



n n 

E llTr([hv, W]p(t)) = £ ll ik) Tr([I n(k)v , W']p'(t)). (55) 

k=l k=l 



Using the fact that 7& = 7^), we can rewrite equation (55) as 

n 

£ t£ (Tr([/ to , W]p(t)) - Tr([/^, W>'(*))) = 0. 
k=i 

Since the coefficients 7^ are all different and the previous equality holds for every I > 0, 
using a Vandermonde determinant type of argument, we obtain 

Tr([I kv , W]p{t)) - Tr([I Ak)v , W']p'{t)) = 0, (56) 

as desired. □ 
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Appendix C: Proofs of the lemmas in Section 5 

Proof of Lemma 5.1 

Proof. 

1. This fact follows easily by applying Lemma 4.4 with W — W — I klVl ,...,k r vri since, if 
Vj 7^ Wj, it holds: 

IklVl,...,k r V r ; IkjWj iIklVl,...,kj{VjWj],...,k r V r - (5^) 

Here we have used the notation [vjWj] — v if [cr Vj ,a Wj ] = ±icr v , and agreed to multiply (57) 
by —1 if the minus sign appears. 

2. By applying Lemma 4.2 to equation (22), we have: 

Tr ([A, I klvl ,..., krVr ] p(t)) = Tr ([A', I klvl ,..., krVr ] p'{t)) . 
Now, we apply Lemma 4.4 to the previous equation to get: 



Tr ( [A, Ik ivi ,...,k r vr\ > hv- k 



Using the Jacobi identity, we have: 

[A, I klVlt ... jkrVr ] , I kv - = A, I kl vi,...,k r V r l Ikv^ 



p(t)) = Tr([[A'J klvl _ krV ]J kv _]p'(t)). 

A, Ikv 



+ 



'-kivi,...,k r v r i 



We have 



kivi,...,k r v T j kv^. 



if k {ki, . . . , k r } 

if 3j with k = kj and v k = Vj 

, 'HklVl,...,kj[vjV%.],...,krVr 

Using the fact that (58) holds for any choice of values v k , and (61) we get: 

Tr {[A, [ik^krvrihv,]} PH)) = Tr {[A', [l klVl ,..., krVr , T kVh \] p'(t)) . 
Thus combining the previous equality with (59) and (60) we get: 



(58) 
(59) 
(60) 

(61) 



Tr ( 



A, I, 



kvr 



i IkiVi,...,k r v r 



(62) 



Tr 



p(t)) = Tr ([[A',1-^-] ,h lVl ,..., krVr } p'(t)) 

Notice that equation (62) is of the same type as equation (58); it is enough to replace A with 
[A, Ifa ] (resp. A' with [A',I kv .]). Thus by applying first Lemma 4.4 and then the Jacobi 
identity we get: 

' p(t))+Tr( 

p'(t)) + Tr ( 
On the other hand, using (61) and (62), we get: 

Tr ([[A,I- kv _], [l klVu ..., krVr ,I TvT ]} p(t)) = Tr ([[A' , I kv _], [l klVl ,..., krVr , I Ivt ]} p'(t)) . 

Thus: 



[A, h 



Tr([[A',I- kVh ], 



-^fcl?Jl,...,fc r ?J r ; Ilvj 



<-kivi,...,k r v r i -LIvt 



IkiVi,...,k r v r 
Ikivi,...,k r v r 



[A, hv- k ], IjvjW p{t) J 

[A',I kv .],I IvT ]]p'(t)). (63) 



Tr ( 



'■klVl,...,krVr 1 



[A,I kv _},I Iv _]]p(t))=Tr( 



'fclDl,...,fc r t>r 



[A', I kv ^\,I] 



Ivi 



p 



which implies (23), as desired. 



□ 
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Proof of Lemma 5.2 

Proof. Both statements are a consequence of Lemma 5.1 (equation (23)). First notice that, 
again by Lemma 5.1, it is enough to prove (25) and (26) for a particular choice of {vj} and 
v. We have, for I > k, 

[Hjz, [Hkx, A}} = -Jkiihzjx- (64) 
I. By applying Lemma 5.1 (equation (23)) to (24) and using (64) we get: 



aTr ( 



kl lI kz,lxl 1 klVl,...,k r V r 



p(t)) = aTr ([-J&hzjxihiv!,...,^] p'(t)) 



(65) 



We may assume, without loss of generality, that k = kj and Vj = x. In this case we have: 

— JkT 



Ikz,lxi Ikivi,...,k r v r 



J—%1 

kl kivi,...,kjy,..., kfi)f * 



Combining the previous equality with (65), equation (25) follows easily. 

2. Using the same procedure, we end up again with equation (65), but now both indices k 
and I are in {hi, . . . , k r }. Assume, for example that k± — k and k 2 = I, and take = fjt 2 = x, 
then equation (26) follows since it holds: 

z,k2Xi Ik\x,k2X,...,k r v r ,k3V3,...,k r V r • 



□ 
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